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Q\ ' Abstract 

o 

, We would comment on the results of the paper "a unified scheme 

for flavored meson and baryons" ( P.C.Vinodkumar , J.N.Panandya, 
• V.M.Bannur, and S.B.Khadkikar Eur. Phys. J. A4(1999)83 ), and point 

out some inconsistencies and mistakes in the work for solving the Dirac 
Qh| equation. In terms of an example for a single particle we investigate the 

reliability of the perturbative method for computing the Coulomb en- 
ergy and discuss the contribution to the wavefunction at origin from the 
Coulomb potential. We conclude that the accuracy of their numerical 
^ ' results needs to be reconsidered. 



In the recent paper , based on a relativistic harmonic confinement model (RHM) with 
a phenomenological parameterization for the residual color electrostatical potential and 
using the residual confined one-gluon-exchange potential (COGEP) for the spin-hyperfine 
interaction, an attempt was made to compose a unified understanding of the spectroscopy 
of hadrons from light flavors to heavy ones. The authors of Ref.[|]] computed the masses of 
mesons and baryons of different flavor combinations for qq,qQ,QQ,qqq,qqQ,qQQ,QQQ 
etc. where q and Q stand for light and heavy flavor quarks respectively, with a unique 
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confinement strength . In this framework they obtained mass differences of pseudoscalar- 
to- vector mesons as well as that for l/2-to-3/2 baryons. They predicted the masses of 
the S-state flavor mesons and baryons which are not measured yet and evaluated leptonic 
decay widths of the vector mesons with no further additional parameters. As shown in the 
tables of Q, their numerical results are in good agreement with the experimental data, 
much of which are even better than that obtained in terms of some successful models. 

The idea is obviously remarkable that one can suppose the effective mechanism for 
both light and heavy quarks to be the same and all hadron states can be well described 
in a unified framework. However, unfortunately, we find that there are some evident 
inconsistencies and mistakes in their work while solving the Dirac equation, using the 
perturbative method to compute the Coulomb energy and omitting the Coulomb potential 
to calculate the wavefunction at origin. The ill-treatment greatly influences the reliability 
of their numerical results. 

1. On the solution of the Dirac equation. 

The basic Dirac equation (3) in which describes the confined single particle state 
of the quark under a mean field potential V(r) is the starting point of their calculations. 
It reads 

[i-fdfi — M q — V{r))^ q {r) = (1) 
where V(r) is a potential of form 

F( r ) = I(l + 7o )AV, (2) 

which includes both scalar and vector parts. 

The wavefunction of quark ip q {r) was expressed as a bispinor form [x<?><A?] which sat- 
isfies coupled equations, 

[E - M q - A 2 r 2 ) Xq = -%a • V0 9 (3) 
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[E + M g ]<f> q = -iff ■ V Xq , (4) 

where M q is the quark mass parameter and the E is the eigenvalue of energy. 
The authors of Ref.(IJ introduced an operator U of the form 

U= - \ 2 ( l. f ^ |. (5) 
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(E+M q )' 2 \ (E+M q ) 



Then ip q can be transformed into a form where the lower component is eliminated 

U% = ( o 9 ) ■ ( 6 ) 
They gave the following normalization 

= ixM = i (7) 

and the upper component Xg satisfies 

[_V 2 + A 2 r 2 (E + M q )] Xq = (E 2 - M 2 ) Xq . (8) 
However, it is very easy to prove that U is not a unitary operator because 

tfU = Utf = ± 1, (9) 
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(e+mY 

unless under the extreme non-relativistic limit. 

Thus, Eq.(^) is incorrect. Instead, we should have 

2 

= mu\\+ {E P +M2) )u\^ q ) 

= fel(i+(^iy)lx 9 ) = i. (10) 

Eq.(|i"o|) gives the normalization condition for Xq- I n fact, we do not need the transforma- 
tion operator U, because E + M q ^ 0, from Eq. (||) we have 



Substituting it into Eq.(||), a corresponding equation similar to Eq.(|8|) for \q can be ob- 



tained immediately. The normalization condition Eq.(lO) is satisfied naturally. However, 
an important difference manifests. The lower component of the ip q is not zero and cannot 
be neglected in the calculation of matrix element, which would have a great influence to 
the numerical results of Ref. |Q] . 

It is not difficult to understand this inconsistency. Instead of U we can take an unitary 



operator U = A /l + j^j^U .i.e. 
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Uo = , 1 a ^ <*™ • (12) 

1 + i^kF V {E+Mq) I 

Using U we can eliminate the lower component of ip q too, but instead of the form of Eq. 

(0) we have 



u *=hwTMtfi°)- <l3> 

Evidently, ip q ip satisfies the following equation: 

U [vfdn -M q - V(r)]U^ p (r) = 0. (14) 

Because the commutation relation between the operator p and r is not zero, the equation 
for tpg P (r) is much more complicated than Eq.(^). 

Actually, it is not difficult to obtain an analytical solution of the Dirac Eq. (|l]) directly. 
According to the general way for solving the Dirac equation with central potential M, let 

<wo - ( % ) . w 

where j is the total angular momentum quantum number, I and I' are the orbital angular 
momentum quantum numbers, £lji m is the well-known spherical spinor and 



Now we insert the expression ( |l6| ) into (JTJ) and can obtain the differential equation for the 
radial parts G and F: 

dG(r) k 



dr r 
dF{r) 



G{r) + [E + M q ]F(r) (17) 



dr r 



F{r) — [E — M q — A r ]G(r), (18) 



where k is a quantum number used frequently in solving Dirac equation with central 
potential. It is defined as 



From Eq. ([L7]) we can obtain 

F = ^ [^ + -g|- (20) 
E + M q l dr r J v ; 

We insert the expression for F into Eq.(|l8|) and obtain a differential equation for G of the 
form 

- + + + M^V<? = (£ 2 - M 9 2 )G. (21) 

It is easy to see that, if we consider k as an angular momentum quantum number I, 
Eq.(21) will be the same as the equation for the reduced radial wave function u(r) = Xq ^ . 



Therefore, we can obtain the solution for E, G and then F by using the similar method 
as that used by the authors of Ref.[jl]]. The proper normalization condition is 

/ (G 2 + F 2 )dr = 1 (22) 
Jo 



2. On the negative energy state 

The potential under consideration is of the form of three-dimensional harmonic oscil- 
lator and we can easily obtain solution for eq. @). The authors of ref.JT) gave the single 
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particle energy as (see Eq. (9) in 

E N = ±y/M* + (2N + 3)n N (q). (23) 

Then they claimed that "following Dirac, the negative energy state is interpreted as an- 
tiparticle". In fact, this statement is incorrect. Because a necessary condition for Eq.(|8|) 
having bound-state- solution is E + M q > and by the definition, 

U N = A(E N + M q )? , (24) 

therefore, in this case Eq,(p3|) only has an unique positive real solution. Taking the 
negative sign in Eq. (|23|), we obtain E = —M q , which is not a solution and moreover 
its absolute value is not equal to the positive solution either. Thus, the negative energy 
solution cannot be interpreted as one for an antiparticle. 

3. On the computation of the Coulomb energy 

In , a residual Coulomb potential 

= (25) 
r 

was introduced and then the Coulomb part of the energy was computed perturbatively 
using the confinement basis. They gave 

e n (qiqj)coui = (N\V coul \N). (26) 

Below we will examine the reliability of the result by this approximate computational 
method. For convenience, we are not going to repeat their calculations, instead consider a 
simpler but reasonable model. Let us put the residual Columbic potential (25) into Eq.(|8|) 
and explore the reliability of the computed results in terms of the perturbative method 
for a single particle state instead of the two-particle-states. A gross estimation for the 
reliability which we are concerning can be obtained accordingly. 



The equation to be solved should have the formtl: 

+ A 2 r 2 ] X n q ew = {E-M q )x n q ew (27) 



V 2 n> e ff 

" l2„2i, new / ttt t\ r \„,neui 



E + M q r ^ q v q ^ q 

where ^(E + M g ) is the dynamical effective mass of the quark. Eq. (|27|) can be rewritten 
as 



-V 2 - - + A 2 r 2 (£ + M g )]x" ew = {E 2 - M 2 q ) X n q ew (28) 



A 

r 

where A = (E + M q )a e Jf '. 

It is obvious that, if we consider yP er = —A a s a perturbative potential, the 0-th 



approximation of Eq. (28) is just the same as Eq.Q. For simplicity, let us take the 
normalized solution Xq f° r Eq.(8) given in Ref. Q as the 0-th order approximate wave 
function and assume Oo = Oi = O. For the 15 and 2S states, the corresponding 0-th 
order approximate eigenvalues are 

emfl =E 2 - M 2 = 30, (29) 

and 

em g = E\ - M 2 = 70, (30) 

respectively. 

We can analytically derive the first-order correction to the eigenvalue {XqW per \Xq) an d 



the results are 

eff 



nn\'- -(/:'„ + ,»,/ In;:-'-'! 2^/") (31) 



and 



em ^ = -(E 1+ m q )a e s ff d\^)- (32) 



v 3 V vr 

Besides, we can also calculate the matrix element of the perturbative potential V p ' 
in the {xq} representation, which reads 



vrr = ( Xq (is)\vnxq(2S)) = ~(e + m^jh^H^)- 



(33) 



In Ref. Q, the a e /f is a running coupling constant depending on both flavor and 
energy state and it is determined by a combination of several complicated relations ( 
Eqs.(13)-(17) of Ref. jlj]). For the ground state of a meson, a gross estimation indicates 
that it is about 0.1 ~ 0.2. For convenience, we take ofj* = 0.15 instead of the running 
coupling constant in our calculation below. 

Using the data given in Q as 

A = 2166 {MeV)l,M u = M d = 82.8MeV, 
M s = S57.5MeV,M c = U28MeV,M b = 4636.6MeF, 

we can obtain the numerical results of em[ j, emfj, and V^ er , we list them below in Table 
1. 

Table 1 



Numerical results of 0-th and 1-st approximate eigenvalues. 



flavor 


Q,{MeV 2 ) 


X(MeV) 


em (0) (MeV 2 ) 


em (1) (MeV 2 ) 


em* 1 '] 


V& er (MeV 2 ) 




u 


46780.8 


69.9696 


140342 


-17076.5 


0.121678 


-6971.44 


0.037256 


s 


65992.9 


139.242 


197979 


-40362.1 


0.203871 


-16477.8 


0.062423 


c 


118144 


446.270 


354432 


-173085 


0.488345 


-70661.5 


0.149524 


b 


209335 


1401.06 


628005 


-723324 


1.15178 


-295296 


0.352659 



As it is well known 0], if the approximation makes sense, the following conditions must 
be respected, 

|em (1) | <C em(°l (34) 

and 

\V?* r \ « emf s - emfl = 40. (35) 

l em (i)i 

It is easy to check from the values in the Table 1 that — }§- is about 0.5 for c quark 

em \s 

and turns to be larger than 1 for b quark. It means that the first condition (|34| ) is 

\v per \ 

broken seriously for c and b quark. Because — njy 2 — tqj is about 0.35 for b quark , the 

em 2s r em is 
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second condition (35) does not hold in the b— case either. The breakdown of the constraint 
conditions would undoubtedly undermine the reliability of the numerical results of Q . 

To understand this problem is not difficult, because the Coulomb potential is a short 
range potential. For the heavy quark it is more important than the confinement potential. 
Our results indicate that it is not appropriate to consider the confinement potential as the 
0-th order and the Coulomb potential as a perturbation for a heavy quark system. 

4. On the wave function at the origin. 

In order to compute the leptonic decay width of vector mesons, the authors of Ref.(l| 
used the radial wave function of meson evaluated at center R^g J (0) (given in eq. (26) of Q). 
It means that they only considered the contribution from the relativistic harmonic mean 
field potential and neglected the residual Coulomb potential. However, in fact the later is 
the potential which makes the main contribution in the short range and must predominate 
the value of the wave function at origin. We would like to give more discussions on this 
aspect by using the simple single-particle example described in Eq. ( |2"g| ) above. The 
square of the wave function at origin for nS states can be obtained by using the well- 
known expression @ 

lx ne Wm 2 = ^_ {nSl dV lnS)i (36) 

where /i is the mass of the particle. For the case corresponding to Eq. fl28|), we have 
2fi = 1. If the Coulomb potential is omitted, the result from Eq. (^) is just the same 
as that computed in Eq.(8) of ref.[[ij. Thus, the corresponding squared wave function at 
origin for IS state is 

lx^(0)| 2 = (x^l|:(^)|4 ) s ) = ( ^)l. (37) 

Including contributions of the Coulomb potential, by using the 0-th wave function of the 
Eq. (|28|) we can approximately estimate the contribution of the residual Coulomb potential 



to the wavefunction at the origin, which reads 

lxS 5 (o)l 2 = (xi°Llx(-r)lxi°L) = A(?). 



(38) 



Still taking a e *f = 0.15 and in terms of the values of Table 1 we can obtain the 
numerical results for |Xg°is(0)| 2 an d Ix^isWI 2 - They are shown in Table 2. 

Table 2 

Numerical results of 0-th and 1-st approximations for the squared wave function at origin for IS state. 



flavor 


lxI, U is(0)| 2 (10 6 A/e^) 


lx^ s (0)| 2 (10 6 Mel/ 3 ) 


u 


1.81709 


1.04190 


s 


3.04454 


2.92494 


c 


7.29277 


16.7825 


b 


17.2004 


93.3574 



The data in Table 2 clearly indicate that the contribution of the so-called residual 
Coulomb potential to the wavefunction is very close to that of the harmonic mean field 
for light flavors, moreover, it turns much larger for heavy quarks b and c. The situation 
is much more serious than for eigenvalues. Therefore, in any case, it is not plausible to 
neglect the Coulomb potential at the 0-th order, even though it is a residual one. 



In summary, we think that the idea to construct a unified scheme for flavored mesons 
and baryons is very desirable and the authors of ref.|l]] have made a new instructive trial. 
However, while solving the basic equation established by the authors and calculating the 
residual Coulomb energy, as well as the wave function at origin there are some inconsistency 
and even mistakes, which affect the reasonability and correctness of their numerical results 
and may destroy the reliability of the solution. Therefore, even though the motivation of 
Ref.[Q| is great, the adopted method for obtaining solution needs to be reconsidered more 
carefully. 
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